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e vertical  temperature  gradient  in  the  tropical  oceans  can  be  used  to  operate  ocean  thermal 
power  plants  (OTPP’s)  that  function  as  heat  engines.  Such  energy  extraction  would  entail  very  large 
scale  pumping  of  both  warm  surface  water  and  cold  bottom  water  through  heat  exchangers,  with 
consequent  disturbance  and  extensive  mixing  of  surface  and  deep  waters.  The  temperature  difference 
between  the  warm  and  cold  inflow  water  should  be  maximized  for  plant  efficiency,  but  entrainment 
and  recirculation  between  the  inflow  and  the  outflow  may  reduce  this  temperature  contrast.  Further, 
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20.  Abstract  (Continued^  ^ Ir  ^ ) 

' horizontal  currents  and  turbulence  carry  the  outflow  water  to  large  distances  and  may  produce  signif- 
icant changes.  These  flows  are  crucial  in  determining  the  thermal  resource  availability  and  the  environ- 
mental impact  of  OTPP  operation.  This  report  provides  background  material  and  brief  assessments  in 
several  areas  of  geophysical  fluid  dynamics  (GFD)  that  bear  directly  on  these  problems.  Relevant 
GFD  research  areas  discussed  include:  turbulence  and  thermal  wakes;  ocean  circulation  and  the 
permanent  thermodine;  air-sea  interaction  and  thermocline  variation;  weather  and  climate  modeling; 
and  marine  ecosystems.  The  report  illustrates  how  each  GFD  area  relates  to  specific  OTPP  problems, 
and  emphasizes  the  multiple  disciplines  in  GFD  that  must  be  considered. 
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GEOniYSICAL  FLUID  DYNAMICS  BACKGROUND 
FOR  OCEAN  THERMAL  POWER  PLANTS 

1.  INTRODUCTION 

We  begin  by  outlining  sone  fundamental  features  of  ocean  thermal 
power  plants  (or  OTPP's).  Next  we  describe  the  problems  of  thermal 
resource  availability  and  environmental  Impact  which  can  be  expected 
In  the  operation  of  such  power  plants.  Then  we  Introduce  the  five  areas 
of  geophysical  fluid  dynamics  (GFD)  that  have  direct  bearing  on  the 
study  of  these  problems.  The  remaining  sections  of  this  report  provide 
partial  surveys  of  the  literature  and  brief  assessments  of  the  research 
methods  in  the  five  GFD  areas. 

Vertical  thermal  gradients  in  the  oceans  can  be  used  to  operate 
heat  engines  of  various  configurations.  In  tropical  waters,  a temperature 
difference  of  about  20°C  exists  between  the  sun-warmed  surface  water  and 
the  deep  cold  water.  The  more  advanced  designs  of  heat  engines  Intended 
to  extract  energy  from  such  a temperature  difference  are  based  on  a closed 
cycle,  employing  a working  fluid  like  ammonia  or  propane.  Heat  is 
supplied  to  the  working  fluid  by  passing  it  through  heat  exchangers 
(evaporators)  in  contact  with  the  warm  surface  water.  This  serves  to 
vaporize  the  working  fluid,  at  a high  pressure;  the  high  pressure  gas 
can  be  used  to  drive  a power  turbine,  provided  it  Is  recondensed  at  low 
pressure  upon  emerging  from  the  turbine.  The  cooling  of  the  gas  leading 
to  condensation  would  occur  in  heat  exchangers  (condensers)  in  contact 
with  cold  water  pumped  up  from  a considerable  depth.  The  working  fluid 
continually  repeats  this  cycle,  effectively  transferring  a considerable 
flux  of  heat  from  the  warm  to  the  cold  water,  while  converting  a small 
Note:  Manuscript  submitted  October  7,  1976. 
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porclon  of  It  to  usable  energy.  With  a 20°C  temperature 
difference,  the  maximum  theoretical  efficiency  of  a heat  engine 

is  about  7 percent,  but  allowing  for  losses  and  for  power  required;  to  drive 
the  pumps,  the  practical  efficiency  would  lie  between  0.5  and  2 percent. 

At  such  low  efficiencies,  ocean  thermal  power  plants  must  pump  vast 
quantities  of  both  warm  and  cold  water  through  the  heat  exchangers  to  produce 
commercially  viable  amounts  of  energy.  At  2 percent  efficiency,  a 500  MWe 
OTPP  must  pump  nearly  2000  cubic  meters  per  second  (or  of  order  1 cubic 
kilometer  per  week)  of  both  warm  and  cold  water  through  its  evaporator  and 
condenser.  Such  a large  scale  movement  of  water  masses  is  characteristic 
of  all  the  OTPP  designs  considered  so  far.  The  OTPP  in  several  design  forms 
can  be  viewed  as  a massive  cold  water  pipe,  300  to  1000  meters  in  length  and 
about  40  meters  in  diameter,  extending  downward  from  the  main  plant  submerged 
near  the  warm  surface.  The  plant,  with  further  water  Inlets  and  outlets, 
accommodates  the  water  pumps,  large  heat  exchanging  evaporators  and  condensers, 
and  electricity  generating  turbine  machinery.  To  supplement  this  brief  summary, 
we  refer  to  the  discussions  of  OTPP’s  by  Lavi  (1973),  Zener  (1973a, 1973b, 1974) , 
and  Heronemus  (1973, 1974a, 1974b) . The  most  recent  design  considerations  and 
engineering  evaluations  of  OTPP's  are  summarized  in  Dugger  (1975),  an  outcome 
of  the  Third  Workshop  on  Ocean  Thermal  Energy  Conversion  (OTEC) . 

Such  power  generation  using  OTPP's  is  environmentally  attractive.  Since 
the  energy  source  is  solar  heating,  it  is  not  readily  depleted.  Nor  is  there 
any  fuel  cost  or  radioactive  danger  or  waste  disposal  problem.  Effects  on  the 
climate  and  on  marine  life  may  perhaps  be  tolerable  or  beneficial,  though  this 
remains  to  be  demonstrated.  The  technical  feasibility  of  OTPP  operation  seems 
close  at  hand. 

The  problems  at  first  sight  are  mostly  in  the  engineering  area,  concerned 
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with  design,  capital  cost,  and  maintenance.  They  Include  the  need  for  cheap 
heat  exchangers,  the  ocean  engineering  problems  of  cold  water  pipes 
and  anchoring,  energy  transmission  from  a floating  OTPP,  and  blofoullng 
of  the  flow  system.  However,  quite  crucial  questions  are  raised  about 
the  thermal  gradient  resource  availability  and  about  the  environmental 
Impact  of  operating  OTPP's.  These  two  areas  require  careful  attention. 

Though  ocean  thermal  power  plants  are  not  accompanied  by  the 
ejection  of  net  heat  into  the  environment,  the  extraction  of  large 
amounts  of  heat  from  localized  and  extended  areas  of  the  upper 
ocean  may  lead  to  a perturbation  In  the  delicate  balance  of  ocean- 
atmosphere  dynamics,  causing  possible  regional  climatic  changes.  A 
group  of  operating  OTPP's  will  erode  the  thermocllne  In  large  areas 
of  the  tropical  ocean,  and  the  knowledge  of  Its  recovery  to  a previous 
or  permanently  altered  state  Is  vital  to  the  assessment  of  resource 
availability  and  environmental  Impact.  To  maintain  efficient  operation, 
OTPP's  must  be  operated  in  a thermocllne  with  a sufficiently  large 
temperature  contrast,  the  exact  value  depending  on  engine  design  and 
power  output.  Hence,  they  must  either  be  anchored  in  a current,  pre- 
sumably close  to  shore  or  in  fairly  shallow  water,  or  be  propelled  at 
some  velocity  to  maintain  availability  of  this  temperature  contrast. 

The  spatial  extent  of  their  "wakes"  or  regions  of  Influence  must  be 
known  in  order  to  minimize  interaction  between  plants.  The  time 
scales  of  the  recovery  of  a thermocllne  must  be  understood  in  order 
to  prevent  the  propelled  or  "grazing"  kind  of  OTPP  from  returning  too 
soon  to  an  ocean  region  where  the  thermal  gradient  is  still  depleted. 


In  this  report  we  will  attempt  to  Identify  the  role  that  geophysical 
fluid  dynamics  (GFD)  plays  In  these  Issues  of  resource  availability 
and  environmental  Impact  of  OTPP  operation.  In  doing  so,  we  will 
provide  brief  summaries  of  the  status  of  research  progress  In  a selection 
of  relevant  GFD  areas.  These  are  certainly  difficult  research  areas, 
since  even  the  simpler  formulations  tend  to  be  nonlinear.  Further, 
progress  in  GFD  research  has  often  resulted  from  a complicated  interaction 
between  observation,  experiment,  analysis  and  numerical  simulation. 

Clearly  also  some  of  the  relevant  GFD  aspects,  such  as  stratified  turbulence, 
are  not  yet  well  understood.  Nevertheless,  the  techniques  of  GFD  must 
be  employed  if  one  Is  to  obtain  realistic  estimates  of  the  effects  of 
large  scale  pumping  by  OTPP’s  or  of  how  natural  variations  can  be 
expected  to  modify  the  efficiency  of  operation. 

The  fundamental  need  Is  for  a calculation  of  the  near-plant  flow, 
which  would  provide  mean  disturbances  of  temperature,  salinity,  and 
biological  nutrients,  together  with  their  turbulent  fluctuations.  The 
calculation  Is  required  first  to  determine  the  mean  temperatures  for 
both  OTPP  Inflows  (warm  and  cold),  as  these  temperatures  may  be  modified 
by  possible  recirculation  of  water  near  the  plant.  Secondly,  this 
caluclation  Is  necessary  to  determine  the  OTPP  impact  from  a far-field 
viewpoint.  In  terms  of  vertical  transport  of  heat,  salt,  and 
biological  nutrients.  Detailed  research  on  OTPP  wakes  has  been 
minimal.  McMlchael  (Zener,  1974,  pp.  8-16)  and  Trimble  (1975, 

Appendix  P)  suggest  that  an  OTPP  can  operate  contlnously  without 
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spoiling  its  ovm  thermal  gradient  resource  even  when  anchored  in  a 
stationary  ocean.  We  suggest  instead  tliat  the  turbulent  outflow 
will  recirculate  locally,  and  thus  ruin  the  resource;  careful  calculations 
are  needed  to  identify  the  maximum  OTPP  pumping  rate  and  the  minimal 
motion  of  an  OTPP  relative  to  its  surroundings  to  avoid  this  recir- 
culation. 

The  turbulent  near-plant  flow  is  determined  by  the  OTPP  inflows 
and  jet  outflows, by  the  plant  design,  and  by  the  ambient  ocean 
temperature,  salinity  and  current  distributions.  This  kind  of  flow  has 
not  been  successfully  predicted  to  date.  But  some  closely  related 
studies  have  been  made  on  thermal  wakes  in  the  presence  of  density 
stratification.  We  therefore  present  in  section  2 a brief  background 
review  of  this  work,  as  specifically  related  to  ocean  thermal  energy 
conversion. 

Studies  of  the  far-field  effect  of  OTPP  operation  on  the  thermal 
gradient  resource  require  an  understanding  of  the  processes  which 
maintain  the  average  vertical  temperature  distribution,  and  its 
seasonal  variation,  in  the  undisturbed  ocean.  Briefly,  the  ocean 
surface  is  cooled  at  the  poles  and  warmed  in  the  tropics,  as  a result 
of  the  balance  of  solar  heating,  infra-red  cooling,  evaporation,  and 
sensible  heat  transfer.  The  cooled  polar  water  sinks,  and  as  a result 
of  deep  circulation  the  ocean  is  filled  with  cold  water,  with  an 
overlaying  warm  buoyant  layer  (except  at  the  poles) . Heat  is  trans- 
ferred poleward  from  this  surface  layer  by  wind-driven  surface  currents, 


and  vertically  downward  (against  the  stratification)  by  turbulent 
mixing.  A weak  general  upwelling  of  cold  water  outside  the  polar  regions 
then  maintains  a statistically  steady  temperature  distribution.  The 
geophysical  background  and  literature  for  these  processes  are  outlined 
in  sections  3 and  4.  There  we  deal  with  ocean  circulation  and  the 
permanent  thermocline,  and  also  with  air-sea  interactions  that  produce 
thermocline  variations.  Clearly  it  is  these  same  processes  which 
determine  the  far-fleld  rate  of  recovery  of  the  ocean  towards  its 
undisturbed  state,  after  the  thermal  stratification  has  been  modified 
by  operation  of  one  or  more  OTPP's, 

Also,  the  operation  of  a large  number  of  OTPP's  in  a limited 
geographical  region  is  likely  to  produce  some  reduction  in  ocean 
surface  temperatures.  This  may  cause  a sufficient  change  in  air 
temperature  and  evaporation  to  produce  a slight  alteration  in  the 
regional  climate.  It  is  Important  to  demonstrate  that  any  such  alteration 
is  within  acceptable  limits.  Section  5 presents  background  material  on 
weather  prediction  methods  using  limited  area  and  general  circulation 
models,  and  discusses  some  of  the  results  related  to  climate. 

It  is  well  known  that  marine  life  thrives  in  upwelling  regions  of 
the  ocean,  and  these  are  among  the  worlds  best  fishing  areas.  This  is 
thought  to  be  due  to  the  vertical  transport  of  biological  nutrients, 
mainly  phosphate  and  nitrate, which  support  plant  life  (phytoplankton) 
and  thus  provide  food  for  zooplankton  and  larger  marine  animals, 

Large  scale  OTPP  operation  will  also  produce  substantial  upward  transport 
of  biological  nutrients,  and  thus  may  be  of  considerable  value  in 
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Assessment  of  the  detailed  biological  impact  of  a particular  j 

K 

region  and  scale  of  OTPP  operation  will  require  a detailed  model  of  *r ! 

the  way  in  which  the  present  distribution  of  phytoplankton,  zooplankton, 

j , 

detritus,  fish  and  dissolved  nutrient  is  maintained.  This  model  I : 

should  include  the  effects  of  horizontal  and  vertical  transport  by  the 

mean  flow  and  turbulence,  the  dependence  of  phytoplankton  growth 

on  light,  and  the  sinking  of  detritus  which  depletes  the  nutrient  in 

the  insolated  surface  layer.  Such  a model  is  discussed  in  section  6, 

which  provides  the  necessary  background  for  determining  the  impact  • 

of  OTPP  operation  on  the  marine  ecology. 


2. 


TURBULENCE  AND  THERMAL  WAKES 


The  flow  in  the  inunediate  vicinity  of  an  OTPP  is  turbulent 
because  of  the  outflow  jets,  the  flow  round  the  plant,  and  the 
buoyant  adjustment  of  the  diffused  outflow.  This  turbulence  is 
three-dimensional,  and  its  study  constitutes  a very  difficult  problem. 
Related  previous  research  has  concerned  two  main  areas,  the  dispersal 
of  heated  effluent  from  a conventional  power  plant  and  the  turbulent 
wakes  behind  submarines. 

In  determining  the  far-field  impact  of  OTPP  operation,  we  need 
to  know  the  vertical  distribution  of  turbulent  transport  properties 
in  the  ambient  ocean.  The  OTPP  induced  turbulence  has  decayed  to 
zero  at  these  large  distances,  but  the  OTPP  Induced  changes  in 
temperature  and  salinity  remain,  and  are  modified  over  long  periods 
by  vertical  turbulent  diffusion  in  the  mixed  layer.  The  relevant 
dif fusivlties  must  be  obtained  from  studies  of  this  ocean  mixed  layer 
or  of  the  similar  atmospheric  boundary  layer  problem.  In  either  of 
these  studies,  the  mean  flow  and  turbulence  properties  vary  mainly 
with  the  vertical  coordinate. 

In  this  section,  we  first  describe  the  different  turbulence 
models,  and  then  discuss  their  application  to  the  three  problems  of 
thermal  effluent  dispersion,  submarine  wakes,  and  planetary  boundary 
layer  modeling. 
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(a)  Turbulence  Models 


The  study  of  turbulence,  and  the  practical  modeling  of  compli- 
cated three-dimensional  random  flows,  is  still  in  an  early  stage  of 
development,  and  constitutes  one  of  the  major  unsolved  problems  of 
fluid  dynamics.  Four  distinct  methods  have  been  used  in  studying 
' these  flows;  of  these  four  we  feel  that  only  two  can  yet  be  applied 

j to  OTPP  flows  at  reasonable  cost  and  with  reasonable  hope  of  partial 

I success. 

j First,  Orszag  (1972  a,b,  1973  a,b)  has  computed  certain  turbulent 

flows  exactly,  using  numerical  representations  of  the  Navier  - Stokes 
I equations,  and  resolving  all  length  and  time-scales  in  the  inertial 

range,  down  to  those  of  the  smallest  eddies,  where  viscosity 

1 dissipates  the  energy.  The  statistical  averages  of  the  resulting  flows 

i 

should  be  independent  of  changes  in  the  initial  conditions  and  in  the 
resolution.  Such  calculations  are  restricted  by  cost  and  computer 

I 

! size  to  relatively  low  Reynolds  numbers  and  to  simple  geometries; 

they  are  a long  way  from  being  practicable  for  OTPP  flows, 
j The  remaining  methods  all  involve  closure  approximations.  The 

‘ second  method  is  the  direct  interaction  approximation  and  the  test  field 

models  of  Kraichnan  (1959,  1971,  1973)  and  of  Herring  and  Kraichnan  (1971). 
j These  have  given  good  results  for  homogeneous  isotropic  problems 

and  even  in  some  more  complicated  situations,  and  have  the  advantage 
I of  not  requiring  a choice  of  free  parameters,  but  they  are  too  costly 

t 
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for  application  to  OTPP  flows. 

The  third  and  fourth  methods  for  calculating  turbulent  flows 

are  phenomenological,  and  both  involve  splitting  the  flow  variables 

into  mean  and  perturbation  parts,  with  only  statistical  average 

products  computed  for  the  perturbation  parts.  The  third  method  is 

» sab~grid  closure,  in  whicli  the  mean  part  of  the  flow  variables  is 

their  smoothed  representation  on  a coarse  mesh,  which  resolves  only 

the  beginning  of  the  inertial  range.  In  the  fourth  method,  statistical  ■ 

closure,  the  mean  flow  is  the  ensemble  average  over  a large  number 

of  flow  realizations,  and  is  therefore  steady  in  many  applications. 

Botii  of  these  closure  methods  include  turbulent  transport  terms  in  the 

mean  flow  equations.  These  terms  are  modeled  in  three  different  ways: 

in  eddy  dlffusivity  models,  the  turbulent  fluxes  of  momentum,  heat,  etc. 

are  determined  from  gradients  of  the  mean  flow  variables,  using  imposed 

eddy  d if fusivit ies  (possibly  functions  of  position).  In  flrst~order 

models  the  dif fusivities  are  calculated  from  the  turbulent  kinetic 

energy  and  a length  scale.  The  turbulent  kinetic  energy  is  determined 

from  a partly  empirical  transport  equation.  In  sub-grid  closure,  the 

! length  scale  is  proportional  to  the  mesh  size;  in  statistical  closure 

I 

I 

models  it  is  either  an  imposed  function  of  position,  or  is  calculated 
algebraically  from  the  mean  flow  variables  and  their  derivatives,  or 

I 
I 

IS  computed  from  an  empirical  transport  equation.  In  second-order  models, 
the  various  turbulent  flux  components  are  each  computed  from  trans- 
; port  equations;  ti»ese  equations  involve  empirical  generation,  decay. 
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and  diffusion  terms. 


Clearly  both  the  sub-grid  and  the  statistical  closure  methods 
are  essentially  empirical,  and  can  become  very  complicated,  with  large 
numbers  of  constants  to  be  determined.  Hopefully  such  constants 
would  have  universal  validity,  but  it  appears  that  changes  are  needed 
for  each  simulation. 

Sub-grid  closure  models  are  discussed  oy  Deardorff  (1967,  1970, 
1974).  These  models  are  always  three-dimensional,  since  the  mean 
flow  variables  must  describe  the  random  turbulent  fluctuations  for 
the  longest  length  scales  in  the  inertial  range.  Depending  on  the 
resolution,  they  can  range  from  being  similar  to  the  fourth  method, 
right  up  to  being  equivalent  to  the  second  method  (with  the  mesh  so 
fine  that  sub-grid  effects  are  negligible).  The  results  can  be 
validated  either  by  agreement  of  the  statistics  with  observation,  or 
by  convergence  of  the  statistical  averages  with  improving  resolution. 

Statistical  closure  models  are  discussed  by  Mellor  and  Vamada 
(1974),  Stuiimiller  (1975),  Launder  (1975),  Biringen  (1975),  Vager 
and  Zilintin  (1968)  and  Lewellen  et  al.  (1573  a,b,  1974  a,b).  The 
statistical  mean  flow  can  be  a function  of  zero,  one,  two  or  three 
space  variables,  ranging  in  complexity  from  homogeneous  turbulence 
to  an  OTPP  flow.  Some  variant  of  the  statistical  models  probably 
offers  the  best  hope  for  a reasonably  economical  and  accurate  solution 
for  the  OTPP  problem.  The  uncertain  question  appears  to  be  whether 
a realistic  turbulent  diffusivity  distribution  ( or  the  equivalent  with 
second-order  closure)  will  be  sufficient  to  suppress  the  instability 
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of  the  mean  flow.  If  stability  cannot  be  insured,  the  more  expensive 
sub-grid  closure  method  will  be  necessary. 

(b)  Thermal  Effluent  Dispersal 

The  dispersal  of  cooling  water  from  a conventional  or  nuclear  power 
plant  has  received  considerable  study.  A 500  Mwe  plant  produces  about 
2 X 10^  cal/sec  of  effluent  heat  flux,  with  a temperature  excess  of 
about  20*^C,  and  a volume  flux  of  10  cubic  meters  per  second  (compared 
with  2000  cubic  meters  per  second  for  our  illustrative  500  Mwe  OTPP) . 

We  are  interested  in  submerged  outlets  (whether  it  be  in  a river,  a 
lake,  or  an  ocean)  since  these  situations  are  the  most  closely  related 
to  OTPP  wakes.  The  influences  of  bottom  topography  and  of  the  shore- 
line should  preferably  be  weak,  since  these  are  unlikely  to  affect 
OTPP  wakes  significantly. 

Fan  (1967),  Koh  & Fan  (1970),  Baumgartner  & Trent  (1970), 

Policastro  & Tokar  (1972),  Jirka,  Abraham  and  Harlemann  (1975),  and 
Dunn,  Policastro  and  Paddock  (1975)  have  given  extensive  literature 
surveys  and  model  studies  concerning  these  problems.  Ahn  & Smith 
(1972)  and  Mangarella  & Van  Dusen  (1973)  have  made  recent  studies  re- 
lated specifically  to  the  oceans.  The  Important  effects,  roughly  in 
Increasing  order  of  time-scale  are:  a)  turbulent  jet  formation,  with 

entrainment  of  cool  surrounding  water;  b)  turbulent  buoyancy  adjustment 
of  the  jet  to  find  its  own  density  level;  c)  laminar  gravitational 
broadening  of  the  thermal  plume;  d)  near-surface  fluctations  with 
tlme-of-day  and  weather;  e)  vertical  and  horizontal  diffusion  of  the 
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thermal  plume  by  ambient  turbulence;  f)  seasonal  surface  changes;  and 
g)  horizontal  advectlon  to  different  climatic  regions.  There  Is  con- 
siderable uncertainty  In  the  theoretical  models,  particularly  with 
regard  to  turbulence,  and  a great  need  for  detailed  comparison  with  data 
from  experiments  and  observations.  Data  compilations  have  been  made 
by  Sllberman  & Stefan  (1970)  and  Tokar  (1971).  Jlrka,  Abraham  and 
Harlemann  (1975)  conclude  on  the  basis  of  comparison  of  observations 
with  model  results  that  the  models  In  their  present  state  of  development 
are  still  unsatisfactory.  The  main  problems  concern  the  turbulence  represen- 
tations, numerical  convergence  of  the  three-dimensional  computer 
models  at  reasonable  cost,  the  complexity  and  consequent  Inflexibility 
of  the  codes,  and  the  wide  range  of  length  and  time-scales  to  be  con- 
sidered. 


(c)  Turbulent  Submarine  Vakes 

A submarine  leaves  a turbulent  wake,  which  entrains  fluid  from 
different  levels  in  a density  stratified  ocean.  The  resulting  collapse 
and  broadening  of  the  mixed  region  generates  internal  waves.  The 
turbulent  wake  Is  closely  related  to  the  jets  and  wake  of  an  OTPP. 

Miles  (1971),  Gran  (1973),  Bell  (1973),  and  Bell  & Dugan  (1973) 
have  considered  the  flow  around  a submarine,  and  the  lee  waves  (body 
waves,  Internal  gravity  waves)  generated,  Chan,  Hlrt  & Young  (1972), 
Lewellen,  Teske  & Donaldson  (1973a,  1974b),  Ko  (1973)  and  Piacsek  & 
Warn-Varnas  (1975)  have  modeled  the  turbulent  entrainment  processes 
and  the  turbulence  decay  behind  the  submarine,  using  increasingly 
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grid  closure  models,  but  this  has  Involved  considerably  greater 
computational  expense.  Orszag  & Pao  (1973)  simulated  the  full 
equations  for  three  dimensional  shear-layer  turbulence,  at  low 
Reynolds  number. 

In  regard  to  the  far-field  impact  of  OTPP  flows,  the  same 
vertical  turbulent  diffusivity  which  influences  the  natural  distri- 
butions of  heat  and  momentum  (and  in  the  ocean  of  salt)  also  influences 
the  disturbance  produced  by  the  plant.  It  may  be  cheaper,  and  of  adequate 
accuracy,  to  calculate  the  vertical  turbulent  diffusivity  distributions 
as  the  ratio  of  the  observed  mean  turbulent  flux  to  the  observed  mean 
vertical  gradient.  These  diffusivity  distributions  can  then  be  used 
to  calculate  the  changes  due  to  OTPP  operation,  eliminating  any  need 
to  compute  dlffuslvlties  corresponding  to  different  weather  conditions 
and  wind  stresses. 
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3.  OCEAN  CIRCULATION  AND  THE  PERMANENT  THERMOCLINE 


The  existence  of  a sufficiently  strong  thermocline  Is  a prere- 
quisite of  successful  OTPP  operation.  The  thermocline  represents  a 
transition  region  In  the  ocean,  where  the  density  Increases  rapidly 
with  depth.  In  general  a density  profile  may  show  several  regions  of 
large  density  gradients.  Including  a diurnal  thermocline  (which 
disappears  at  night),  a seasonal  thermocline,  and  a permanent  one. 

Ater  large  storms  new  gradient  regions  may  appear,  and  old  ones  dis- 
appear or  get  displaced.  The  maintenance  of  tne  permanent  thermocline 
depends  on  both  the  ocean  currents  and  turbulent  entrainment  due  to 
the  global,  climatic  wind  systems.  The  time-varying  components,  on  the 
other  hand,  are  mostly  determined  by  the  atmospheric  weather  patterns 
and  solar  heating  through  the  air-sea  Interaction  processes.  In  section 
3 we  will  discuss  the  ocean  currents  and  the  permanent  thermocline. 

In  section  4,  the  variable  components  and  the  air-sea  Interaction. 

In  each  case,  we  will  give  a brief  Introduction  to  the  present  under- 
standing of  the  dynamics,  some  of  the  more  successful  modeling  efforts, 
and  some  observational  highlights. 


(a)  Thermocline  Characteristics 


The  permanent  or  main  thermocline  represents  a transition  between 
the  warm  surface  water  and  the  cold  "abyssal  water"  of  the  ocean.  The 
deep  layers  of  the  ocean  are  filled  with  very  cold  water,  which  sank 
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at  the  poles,  and  is  rising  very  slowly  over  the  rest  of  the  ocean.  The 
upper  layers  are  warmed  at  the  surface  by  solar  heating  and  heat  conduction 
from  the  atmosphere,  and  this  heat  diffuses  downward  by  wind  generated 
turbulence. 

This  upwelling  has  an  estimated  global  average  vertical  velocity  of 
2 cm/day  or  7 meters  per  year  (Wyrtkl,  1961)  which  just  balances  the 
average  downward  entrainment  by  storm-driven  eddy  diffusion.  The  seasonal 
and  spatial  variations  in  thermocline  thickness  and  depth  reflect  the 
variations  of  wind  and  surface  temperature,  and  their  effect  on  the  eddy 
diffusivity  distribution  and  the  horizontal  fluxes  in  the  upper  layer, 
in  a way  that  is  not  yet  fully  understood.  The  latitude  variation  of 
thermocline  depth  in  the  Atlantic  Ocean  is  shown  In  Figure  1;  the 
shallow  thermocline  at  the  equator  is  due  to  the  combined  effects  of 
strong  upwelling  and  low  turbulent  diffusivity. 


E 

I 

Fig.  1 — Variation  of  the  depth  of  the  thermocline  in  the  Atlantic  Ocean 
with  latitude  (mean  values  for  the  f>ntire  ocean)  [Defant,  1961] 
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In  Figure  2 are  shown  average  temperature  profiles  for  February 
and  August  In  the  Gulf  of  Mexico.  These  representative  tropical  profiles 
show  a seasonal  variation  above  350  meters  depth.  The  temperature  profile 
below  this  depth  Is  reasonably  Invariant  to  the  season,  and  Is  associated 
with  the  permanent  thermocllne.  In  regions  of  the  ocean  where  strong 
warm  or  cold  currents  flow,  the  main  thermocllne  will  exhibit  a much 
sharper  spatial  variation;  In  some  cases  It  will  coincide  with  the  lower 
or  upper  boundaries  of  the  current.  Figure  3 represents  temperature 
cross  sections  In  the  Gulf  Stream,  again  In  the  months  of  February  and 
August;  the  rapid  spatial  variation  of  the  temperature  between  500  and 
1200m  Is  clearly  visible. 


(b)  Ocean  Current  Systems 

The  most  important  factors  in  determining  the  circulation  of  the 
world's  oceans  are  the  predominant  wind  patterns,  the  pole  to  equator 
temperature  difference,  and  the  topography  of  the  continents  and  the 
ocean  bottom  (Defant,  1961;  Von  Arx,  1962;  Phillips,  1966;  Neumann  & 

Pierson,  1966;  Kraus,  1972).  In  this  sense,  ocean  movements  differ  from 
the  winds,  which  are  constrained  much  less  by  the  mountains  and  are 
driven  by  temperature  differences  to  a dominant  extent. 

One  of  the  challenges  of  modern  oceanography  Is  to  separate  the 
wind  driven  part  of  the  currents  from  the  density  driven  part.  The  depth 
to  which  the  effects  of  the  wind  penetrates  Is  limited,  and  except  for 
topography,  could  not  couple  to  the  deeper  layers.  As  it  is,  however,  wind- 
driven  surface  fluxes  Interact  with  the  continents  and  lead  to  well-known 


18 


upwelllng  areas.  Some  of  the  deeper  cold  currents  of  the  ocean  are 
simply  the  result  of  mass  conservation,  returlng  the  warm  water  blown  by 
the  wind  in  the  other  direction. 

Because  of  topography  the  currents  in  each  of  the  world's  ocean 
basins  tend  to  form  a self-contained  system.  The  North  and  South  Atlantic, 
the  North  and  South  Pacific,  the  Indian  Ocean,  and  the  two  artic  regions 
all  have  well  established  current  patterns.  The  current  system  of  the 
North  Atlantic  and  the  Gulf  of  Mexico  are  given  in  Figure  4,  showing  also 
the  sea-surface  temperatures  and  how  well  they  relate  to  the  currents. 

In  general.  In  each  ocean  basin  the  surface  circulation  is  poleward 
to  the  west.  In  the  North  Atlantic  and  the  North  Pacific  the  gyres  are 
clockwise,  with  the  westward  North  Equatorial  Current  leading  into  the 
poleward  Gulf  and  Kuroshio  Currents  respectively.  This  picture  is  con- 
sistent with  the  general  atmospheric  wind  patterns,  with  the  trade  winds 
driving  the  westward  Equatorial  Current,  and  the  mid-latitude  westerlies 
driving  the  Gulf  and  Kuroshio  eastward,  and  establishing  the  essentially 
wind-driven  character  of  the  water  movements  above  the  thermocline. 

Stommel  (1965)  has  discussed  the  Gulf  Stream  in  great  detail. 

A special  comment  must  be  made  about  the  currents  near  the  equator. 
The  so-called  North  and  South  Equatorial  Currents  are  actually  displaced 
north  of  the  equator  such  that  the  South  Current  is  on  the  equator  and 
the  North  Current  is  on  15°N.  In  between  flows  the  Equatorial 
Countercurrent  in  an  easterly  direction;  this  current  has  a high  velocity 
core  called  the  Cromwell  Current. 
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FICUKE  vii-2.  The  averof’e  surface  temperature  (’C)  of  the  North  Atlantic 
Ocean  in  April.  Compare  the  position  of  the  isotherms  in  this  figure  with 
the  current  paths  depicted  in  Figure  VI 1-3. 


The  deep  circulation  of  the  ocean  consists  of  relatively  strong 
jet-like  currents  that  move  along  the  continental  boundaries,  carrying 
cold  arctic  water  towards  the  equator.  Among  these  are  the  Labrador 
Current  and  the  Brazil  Current.  The  overall  picture  is  rather  complicated 
however,  as  Figure  5 demonstrates  in  showing  how  the  polar  cold  water 
moves  poleward  beneath  the  main  thermocline. 

The  currents  that  are  of  particular  interest  to  the  OTEC  effort 
are  the  Gulf  Stream,  the  North  Equatorial  Current  and  the  undercurrents, 
since  these  will  be  involved  in  the  resource  assessment  and  environmental 
impact  calculations  for  OTPP  operations.  The  sea-surface  temperatures  of 
the  western  North  Atlantic  are  shown  in  Figure  6 for  both  a typical  month 
of  February  and  September.  The  resource  availability  becomes  clearer  in 
Figure  3 in  the  vertical  slices  taken  through  the  Gulf  Stream  for  the 
same  two  months;  these  displays  of  temperature  with  depth  and  cross- 
stream coordinate  readily  shown  the  temperature  contrast  across  the 
thermocline. 

Finally,  we  show  in  Figure  7 the  run  of  horizontal  velocity  in  a 
similar  slice  across  the  Gulf  Stream;  as  is  true  also  of  the  Kuroshio, 
the  currents  in  the  deep  are  much  stronger  than  expected  from  simple 
wind-drift  or  convective  circulation  theories. 

The  Initial  theoretical  studies  of  ocean  circulation  focused 
on  wind-driven  motions  of  a homogeneous  ocean,  and  these  we  will  not 
review  here.  The  more  recent  modeling  studies  fall  into  three  general 
categories:  a(  simulations  of  the  global  ocean  circulation;  b)  computations 
for  single  ocean  basins;  c)  treatment  of  special  regions  such  as  upwelllngs 
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Fig.  5 — The  fairly  intricate  deep  sea  circulation  of  cold  water  and  the  wanner  surface 
currents  of  the  Atlantic  Ocean  are  shown  [Defant,  1961] 


and  currents. 

Global  simulations  of  a baroclinic  ocean  have  been  undertaken  by 
Crowley  (1970a, b)  and  Bryan  (1969).  Although  these  numerical  models  employed 
rather  coarse  computational  grids  in  the  horizontal,  they  did  succeed  in 
replicating  many  of  the  large-scale  observed  features  of  the  ocean  cir- 
culations. Those  features  include  the  western  boundary  currents,  like 
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Fig.  6 — Contours  of  sea-surface  temperature  in  the  western  North  Atlantic 
are  shown  for  February  and  August  (Stommei,  1965] 


Fig.  7 — Velocity  profile  across  the  Gulf  Stream  (off  Chesapeake  Bay, 
20-22  April  1932).  Values  displayed  in  cm/sec.  [Defant,  1961] 


the  Kuroshio  and  Gulf  Stream,  and  the  Equatorial  Undercurrent.  The 
intensities  of  each  of  these  main  currents  in  the  models  were  not  too 
different  from  the  observed  values,  but  this  resulted  partly  from  tuning 
of  the  models.  Simulation  of  single  ocean  basins  have  been  made  for  the 
Arctic  Ocean  by  Galt  (1973)  and  Semtner  (1973);  for  the  Pacific  Ocean  by 
O'Brien  (1969),  Alexander  (1973)  and  Huang  (1974);  and  for  the  Atlantic 
Ocean  by  Friedrichs  (1972). 

The  spatial  resolution  of  the  single  basin  models  was  several 
times  better  than  that  of  the  global  models,  and  they  provide  slightly 
better  descriptions  of  the  main  currents. 

Recent  modeling  efforts  have  concentrated  on  understanding  the 
newly  discovered  mesoscale  ocean  eddies  (of  ~ 300  km  diameter).  But 
these  studies  mainly  aim  to  understand  the  baroclinic  instability  and 


(b)  Seasonal  Tbermocline 


Figure  9 snows  tlie  seasonal  variation  of  two  representative  temp- 
erature profiles  in  mid-latitudes.  A shallow  seasonal  thermocllne  forms 
in  June  in  both  cases,  and  deepens  and  weakens  until  It  disappears  in 
February,  effectively  merging  into  the  main  thermocllne  below.  The 
adjoining  panel  shows  the  monthly  average  of  sea-surface  temperature  and 
downward  neat  flux  at  a particular  station.  Figure  2 shows  the  seasonal 
variation  in  the  subtropics  between  the  months  of  February  and  August.  Tlie. 
disappearance  of  the  strongly  mixed  layer  (down  to  300  m in  February)  in 
August  is  striking. 

The  difficult  problem  of  simulating  the  seasonal  thermocllne  was 
attempted  by  Kraus  & Turner  (1967)  and  by  Turner  (1969).  They  constructed 
a simple  theory  in  which  all  the  kinetic  energy  of  wind  stirring  was  assumed 
to  change  tne  potential  energy  of  the  system,  and  using  a saw-tooth  form  of 
the  heating  function,  calculated  temperature  and  depth  of  the  mixed  layer 
as  a function  of  time.  Tney  argued  that  in  reality  all  effects  of  stirring 
put  in  near  the  surface  propagate  down  without  being  affected  by  rotation 
or  advectlon. 


(c)  Thermocllne  Erosion  During  Storms 


The  strong  erosion  of  tliermoclines  during  the  passage  of  hurricanes 
has  been  tne  subject  of  several  studies,  among  these  the  analytical  studies 
of  Geisler  (1970),  the  numerical  studies  of  O'Brien  & Reid  (1967),  O'Brien 
(1967),  and  Madala  & Piacsek  (1974),  and  the  ocean  observations  of  Leipper 
(1967)  (see  Figure  12  for  some  of  his  measurements  before  and  after  hurricane 
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Fig.  8 (left  panel) 


FlO.  J.  Santple  KHindini;  data  for  the  iho  ut%  of  Hau.  On  tlic  left  are  the  actual  r(:) 
curvet  while  ihoic  on  the  hi'ht  are  diiTerencet  from  the  iem|)erature  at  M.5  m. 


Fig.  8 ~ Daytime  heating  cautes  the  formation  of  a diumai  thermodine  in  these  time 
histories  of  water  temperature  with  depth,  shown  at  selected  points  in  the  left  panel  and 
as  profiies  in  the  right  panels  [Kaiser  and  Williams,  1974] 


Fig.  9 — In  the  left  panels,  the  seasonal  variation  of  temperature  with  depth  shows  the  buildup  of  a 
seasonal  thermocline  during  the  summer  months  and  its  decay  during  the  winter.  The  right  panel 
displays  sea-surface  temperature  and  downward  heat  flux  throughout  a year  at  a given  location. 
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passage).  Both  upwelllng  due  to  a diverging  wind-driven  flux  in  the  surface 
layer  and  mixing  due  to  turbulent  entrainment  are  shown  to  contribute  to 
the  development  of  the  new  temperature  structure.  The  recovery  of  this 
tliermocline  to  a shape  resembling  pre-storm  conditions  is  shown  by  the 
observations  to  take  approximately  two  months;  the  models  did  not  simulate 
tne  recovery  phase.  A cyclonic  current  system  is  also  set  up  after  the 
hurricane  has  passed  by  the  large-scale  perturbation  of  the  density  field; 
it  disappears  in  a few  weeks. 

The  total  heat  carried  from  the  warm  surface  water  to  the  colder 
deep  water  because  of  the  hurricane  is  approximately  A500  cal/cm^,  spread 
over  an  area  of  1.5  x lO"*  kra*^.  The  sea-surface  temperature  is  lowered  by 
5°C  or  more.  The  numerical  models  successfully  simulated  these  observations 
both  the  axisymmetric  model  of  O'Brien  L Reid  and  the  three-dimensional 
mooel  of  Madala  & Piacsek;  the  latter  included  the  movement  of  the  hurricane 
and  the  results  showed  a more  realistic  spatial  and  long-time  behavior  of 
the  temperature  field. 


(d)  One-dimensional  Models  of  Wind  Mixing 

Denman  & Miyake  (1973)  and  Denman  (1973)  have  investigated  both 
observationally  and  numerically  the  rate  and  extent  of  the  wind-induced 
deepening  of  the  mixed  layer  during  the  passage  of  several  weather  dis- 
turbances in  the  mid-latitudes.  They  also  observed  the  formation  of  a 
shallow  layer  of  warm  water  under  conditions  of  low  winds  and  intense 
solar  heating.  Their  one-dimensional  numerical  model  successfully  simu- 
lated the  observed  behavior  during  a twelve  day  period,  despite  neglecting 
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horizontal  and  vertical  advection  in  the  ocean.  They  concluded  that  for 
their  case  of  spring  development  of  the  summer  thermocline,  these  processes 
were  not  important.  Pollard  et  al . (1973)  have  constructed  an  integral 
model  to  describe  the  response  of  the  upper  ocean  to  an  imposed  stress 
and  to  heating  or  cooling.  They,  too,  found  that  on  time  scales  of  a day 
or  less,  stress  dominates  over  heating  and  cooling  in  producing  eddy 
diffusivity. 

Mellor  & Durbin  (1975)  constructed  a model  for  the  mixed  layer 
and  the  thermocline  using  mean  flow  equations  and  the  turbulent  Reynolds 
stress  equations.  They  obtained  excellent  agreement  with  observations 
provided  by  Denman  & MiyaRe. 

(e)  Air-sea  Interaction 

As  a whole,  the  oceans  receive  as  much  heat  as  they  give  up; 
however  the  heat  gains  and  losses  have  a latitudinal  distribution  from 
equator  to  pole.  More  than  half  of  the  excess  heat  absorbed  by  the  tropical 
oceans  is  given  to  the  atmosphere  tiirough  evaporation,  and  the  atmosphere 
transports  the  heat  (as  latent  heat  of  evaporation)  toward  polar  regions. 

The  large  ocean  currents  also  contribute  to  tue  poleward  flow  of  heat;  the 
best  estimate  of  their  share  is  about  30  percent.  The  latitude  variation 
of  tne  various  forms  of  heat  transport  are  summarized  in  figure  10. 

Recent  studies  (Simpson,  1969) have  indicated  that  in  mid-latitudes 
the  direct  effects  of  ocean  heating  upon  cyclone  growth  is  small,  but  the 
indirect  effect,  via  condensation  as  rainfall  and  corresponding  latent  heat 
release,  is  large.  As  we  shall  see  further  in  section  5.b«  sea-surface 
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Fig.  10  — Averages  and  bounds  of  various  estimated  heat 
transports  poleward  [Neumann  and  Pierson,  1966] 


temperatures  associated  with  ocean  currents,  or  with  areas  of  temperature 
anomalies  in  the  ocean,  have  a profound  effect  on  regional  climates  in  the 
respective  areas  of  the  world.  The  daily  incidence  of  solar  radiation,  and 
its  seasonal  variation,  together  with  the  surface  wind  stress  and  wave- 
generated  turbulence,  essentially  determine  the  location  and  strength  of 
the  diurnal  and  seasonal  thermoclines. 

An  intimate  knowledge  of  the  air-sea  exchange  processes,  in 
particular  the  transmission  of  heat  and  momentum  and  evaporation  rate,  is 
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necessary  if  we  are  to  understand  and  successfully  model  tiiemiocline 
behavior  and  climate  impact  of  sea-surface  temperature  changes. 

The  motions  and  their  accompanying  transports  in  the  air-sea 
boundary  layer  divide  into  two  groups,  depending  on  whether  the  water  is 
warmer  or  colder  than  the  air.  The  first  situation  generally  occurs  at 
night  and  in  the  fall.  Strong  evaporation  and  infrared  back-radiation 
lead  to  surface  cooling  that  in  turn  Jeads  to  convective  overturning  in 
the  ocean.  The  convection  enhances  the  tlow  of  heat  from  the  lower 
layers  to  the  surface,  at  the  same  time  mixing  up  the  water  and  removing 
any  shallow  therraociine.  The  flow  of  lieat  upwards  makes  the  overlying 
air  boundary  layer  unstable  by  reducing  its  stability.  So,  both  an 
increase  in  wind  shear  turbulence  and  in  buoyancy  driven  convective 
motions  will  result  (Delnore,  1972;  Raiser  & Williams,  1974;  Ostapoff  & 
Worthem,  1974). 

When  the  air  is  warmer  than  tiie  water,  the  situation  is  very 
stable  in  both  the  water  and  the  air,  and  results  in  accumulation  of 
warm  surface  waters,  aldii.g  the  diurnal  and  seasonal  thermocline  formation. 
Evaporation  is  a function  of  air  and  sea  temperatures  and  of  the  relative 
iiumidity  of  the  air,  and  of  wind  speed.  Observations  indicate  that  the 
maximum  rate  of  evaporation  in  the  oceans  occurs  between  10°  and  20° 
latitude,  this  being  because  the  cloud  cover  is  usually  greater  and  the  wind 
speed  smaller  at  the  equator.  Bowen's  ratio,  R,  is  defined  as  the  ratio 
of  sensible  to  evaporation  heat  loss  from  the  ocean  to  the  atmosphere. 
Extensive  field  measurements  in  various  oceans  established  that  in  turbulent 
air-stia  interfaces  R has  a value  which  is  of  order  0.1;  its  value  is 
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practically  independent  of  wind  speed,  and  dependent  only  on  the  ocean 
temperature  and  the  air  temperature  and  humidity. 


We  have  concentrated  our  attention  here  on  heat  flow  processes, 
since  these  are  most  crucial  to  OTPP  operation.  A large  body  of  literature 
exists  on  the  generation  and  maintenance  of  surface  waves,  and  their  associ- 
ated turbulence  and  stress;  these  will  not  be  detailed  here.  Suffice  it 
to  say  that  turbulence  in  both  the  air  and  the  water  can  be  nandled  reasonably 
adequately  by  first  or  second-order  closure  metnods  (Deardorff,  1970; 

Lewellen  et  ai,  1974a;  Mellor  6<  Durbin,  1975)  or  by  analytical  prescriptions 
(Munk  & Anderson,  1948;  Monin  & Obukhov,  1954;  Pandolfo,  1966;  Yoshihara, 

1968).  An  extensive  discussion  of  the  many  aspects  of  air-sea  inter- 
actions is  provided  by  Kraus  (1972). 
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i 5.  WEATHER  AND  CLIMATE  MODELING 

Tneoretical  models  concerned  with  the  large-scale  dynamics  of  tlie 
atmosphere  can  be  discussed  in  four  broad  groups.  Although  all  these  models 
respond  in  some  measure  to  the  ocean  surface  temperature,  little  progress  has 
been  made  in  predicting  atmosphere  and  ocean  dynamics  simultaneously.  AJl 
these  models  play  some  role  if  one  wishes  to  estimate  what  effects,  if  any, 
multiple  OTPP  operation  may  have  on  weather  and  climate.  The  first  three 
models  deal  with  scales  smaller  than  global,  and  these  regional  treatments 
are:  a)  the  limited  area  forecast  models  (or  LAi'l's)  for  the  atmosphere; 
b)  studies  correlating  sea-surface  temperature  anomalies  with  regional 
atmosphere  effects;  c)  hurricane  simulation  models.  Finally  in  d)  we 
briefly  outline  the  global  weather  models  or  general  circulation  models 
(GCM's),  and  discuss  what  may  be  inferred  from  tnese  models  about  climate 
prediction. 

(a)  Limited  Area  Forecast  Models 

The  group  operation  of  100  to  1000  OTPP's  calls  for  a region  of 
typical  diameter  1000  to  3000  km.  The  resulting  sea-surface  temperature 
and  heat  flux  anomaly  will  affect  the  atmospheric  motions  above  this  re- 
gion, and  may  produce  some  climate  changes  both  in  the  operational  region 
and  at  greater  distances. 

Wiien  studying  atmospheric  dynamics  over  a small  portion  of  the 
globe,  limited  area  models  are  embedded  in  a large-scale  (usually  global)  flow. 
lam's  employ  a very  fine  mesh,  with  grid  size  10  to  50  times  smaller  than  the  grid  size 
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associated  with  global  circulation  models.  Such  good  spatial  resolution 
is  necessary  if  local  weather  features  are  to  be  adequately  modeled.  The 
boundary  conditions  for  the  LAM  are  taken  to  be  values  determined  either 
from  a GCM  or  from  field  observations,  and  have  prescribed  time  variation. 
Usually  the  large-scale  flow  is  assumed  to  be  unaffected  by  the  LAM  flow. 

Birchfield  (i960),  Wang  & Halpern  (1970),  Shapiro  & O'Brien  (1970), 
and  Asselin  (1972)  all  integrated  barotropic  vorticity  or  primitive  equation 
lam's,  and  Hill  (1968)  a two-level  baroclinic  version.  Considerably  more 
elaborate  multi-layer  LAM's  are  those  of  Williamson  6 Browning  (1974)  at 
the  National  Center  for  Atmospheric  Researcli  (NCAR) , and  of  Hoverraale  (1974) 
at  the  National  Meteorlogical  Center  (NMC). 

Hill,  Wang,  & Halpern  obtained  boundary  conditions  for  the  LAM's  by 
interpolation  of  the  data  produced  by  a coarse  GCM  integration;  this  in 
time  produced  unwanted  small-scale  spatial  oscillations  in  the  solutions 
at  the  outflow  edges  of  the  LAM  regions.  Shapiro  & O'Brien  and  Asselin 
applied  similar  interpolated  boundary  conditions  at  inflow  points,  but 
advected  quantities  out  of  the  region  with  a quasl-Lagrangian  scheme  in 
order  to  suppress  such  oscillations.  Oliger  & Sundstrohm  (1975)  suggested 
a suitable  viscous  term  to  suppress  the  spatial  oscillations  in  the  LAM, 
and  Williamson  & Browning  and  Madala  (1973)  have  used  it  successfully.  In 
summary,  these  various  versions  of  LAM's  iiave  tended  to  demonstrate  a good 
potential  for  forecasts  of  a day  or  two. 
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(b)  Sea-surface  Temperature  Anomaly  Effects 

First  we  shall  consider  the  observational  evidence  for  the  per- 
manent and  temporary  effects  on  the  climate  by  sea-surface  temperature. 

Next,  we  will  examine  the  previous  efforts  to  computationally  simulate 
these  anomalies  phenomena. 

The  crucial  role  of  the  oceans  in  climate  formation  has  long  been 
recognized.  The  rather  mild  climate  of  northwestern  Europe  is  attributed 
to  the  nearby  passage  of  the  Gulf  Stream;  the  rather  cool  and  dry  climate 
of  California  in  the  summer  is  attributed  to  strong  upwellings  of  cold 
water  off  the  Pacific  coast.  More  recently,  scientists  have  also  observed 
and  tried  to  model  associations  between  sea-surface  temperature  anomalies 
and  changes  in  long-range  weather  patterns.  Bjerknes  (1966,  1969)  has 
suggested  that  anomalies  in  the  tropical  Pacific  influence  pressure  and 
circulation  patterns  in  the  higher  latitudes.  Namias  (1970,  1973)  has 
observed  seasonal  persistence  of  a large  (4000  km  by  2000  km)  sea-surface 
temperature  anomaly  of  +6*^0  in  the  North  Pacific,  and  of  its  effects  on 
the  atmosphere.  Ratcliffe  & Murray  (1970)  have  made  similar  observations 
on  correlations  between  large  (2000  km  by  1000  km)  North  Atlantic  sea-surface 
temperature  anomalies  of  approximately  -2°C  south  of  Newfoundland  and 
atmospheric  circulation  patterns  of  western  Europe,  extending  over  the 
time  period  1888  to  1970,  (see  figure  11). 

In  general,  the  effects  of  such  anomalies  are  felt  downwind, 
carried  by  the  global  wind  systems  (easterly  trade  winds  in  the  tropics, 
westerlies  in  the  mid  latitudes.  Spar  (1973  a,  b,  c)  has  simulated  the 
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^■•^!urc  1.  (a)  Mean  pattern  of  sea  surface  temperature  anomaly  (deg  C)  in  September  based  on  years  with 

* *>gnificant  negative  SST  anomaly  area  off  Newfoundland.  Positions  of  anomaly  centres  and  last  two  figures 

of  year  shown  : 

■ negative  anomaly  > 2 deg  C 
A negative  anomaly  between  1 and  2 deg  C 
X position  and  intensity  less  certain. 

■b)  Kfean  surface  pressure  anomaly  (mb)  in  October  for  the  years  with  a negative  SST  anomaly  centre  off 
Ntfwfoundbnd  in  September  as  shown  in  Fig.  I (a).  Stippling  shows  area  where  mean  anomaly  is  significant 

at  5 per  cent  level  according  to  ‘ J ' test. 

Note  : in  this  and  similar  Figures  rcro  anomalies  arc  shown  as  broken  lines. 

Fig.  11  — Sea-surface  temperature  anomaly  off  Newfoundland  is  shown  in  upper  panel, 
which  may  influence  atmospheric  circulation  pattern  over  western  Europe  [Ratcliffe  and 
Murray,  1970| 
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the  monthly  and  seasonal  effects  of  the  North  Pacific  warm  anomaly  re- 
ported by  Namias  (1973)  using  the  Hintz-Arakawa  two-layer  GCM.  In  ad- 
dition to  investigating  the  observed  situation,  he  placed  the  anomaly 
in  the  Southern  Hemisphere.  His  experiments  showed  significant  inter- 
hemispheric  effects  after  about  one  month,  phase  shifts  of  one  to  two 
weeks  in  major  cyclone  development,  stronger  reactions  to  anomalies  in  the 
winter  than  in  the  summer,  and  marked  difference  in  behavior  when  the 
anomaly  was  placed  in  Che  Southern  Hemisphere.  Kowntree  (1972)  has  used 
the  nine-level  hemispheric  model  developed  at  the  Geophysical  Fluid 
Dynamics  Laboratory  of  NOAA  for  anomalies  in  Che  tropical  Pacific,  and 
obtained  results  similar  to  chose  of  Spar.  Houghton  et  al  (1974)  have 
simulated  the  Newfoundland  anomaly  using  the  NCAR  six-level  GCM.  They 
found  good  agreement  with  the  observed  data  in  the  case  of  the  warm 
anomaly;  Che  response  in  the  case  of  the  cold  anomaly  was  not  as  satisfactory. 

There  do  not  appear  to  be  studies  using  limited  area  forecast  in 
connection  with  sea-surface  temperature  anomalies.  This  must  be  due  partly 
to  the  fact  that  the  respective  anomalies  were  far  removed  from  any 
population  centers  or  even  land  areas,  so  there  was  no  interest  to  explore 
their  effect  on  neighboring  regions. 

The  above  observations  suggest  that  anomalies  of  magnitude  2 to 
6°C  and  covering  areas  ^ 10^  km^  may  cause  global  weather  changes.  That 
sea-surface  temperature  anomalies  affect  neighboring  areas  as  well  is 
clear  from  the  climates  of  western  Europe  and  California,  if  we  regard  Che 
Gulf  Stream  and  the  east  Pacific  upwelling  as  "permanent  anomalies". 


associated  with  a passing  hurricane  could  become  a serious  impediment  to 
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OTPP  operations  (see  figure  12).  The  hurricane  models  that  have  been 
employed  for  research  or  operational  forecast  have  all  been  confined  to 
atmospheric  motions  only,  despite  a strong  energy  and  momentum  exchange 
with  the  ocean. 

Some  of  the  present  atmospheric  models  are  those  of  Anthes  (1972), 
Mathur  (1974),  Madala  (1973),  Madala  & Piacsek  (1975),  and  Kurihara  & Tuleya 
(1974).  The  Important  features  with  which  one  can  distinquish  these  models 
are:  a)  the  type  of  mesh  used;  b)  the  number  of  layers  in  the  vertical; 

c)  the  manner  in  which  latent  heat  released  by  cumulus  clouds  is  parameterized; 

d)  air-sea  interaction  and  boundary  layer  treatment;  e)  size  of  domain 
treated,  including  the  use  of  an  f-plane  (constant  Coriolis  force)  or  a 
3-plane  (North-South  variation  of  Coriolis  force);  f)  Imposed  large-scale 
flow  and  boundary  treatment  of  the  fine  mesh. 

The  four  models  give  results  that  differ  mainly  in  the  finer 
details.  A general  statement  that  can  be  made  of  all  four  is  that  they 
successfully  simulated  the  magnitude  of  the  central  pressure  minimum,  the 
rate  of  Intensification  of  the  storm,  the  spiral  cloud  bands  and  associated 
confluence  and  diffluence  lines  in  the  stream  function,  the  organized 
bands  of  clouds  near  the  eye,  the  downward  motion  in  the  eye,  and  the 
cyclonic  outflow  in  the  upper  troposhere.  Mathur 's  and  Madala 's  models 
are  on  a 3-plsne,  allowing  the  hurricane  to  move,  and  all  except  Kurihara 
have  three  layers  in  the  vertical,  while  he  used  eleven.  All  have  used 
explicit  time  differencing  techniques,  except  Madala,  who  used  a semi- 
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implicit  scheme  to  overcome  fast  gravity  waves.  Also,  Mathur  used  two, 
and  Madala  three,  nested  grids,  the  finest  mesh  located  over  the  central 
part  of  the  hurricane.  In  addition,  Madala's  finest  grid  was  movable, 
following  the  path  of  the  hurricane,  Mathur 's  and  Kurihara's  were  fixed 
on  the  eye  of  the  stationary  hurricane. 

(d)  Global  Weather  and  Climate  Modeling 

We  have  seen  in  section  (b)  that  sea-surface  temperature  anomalies 
with  a modest  dimension  of  a few  thousand  kilometers,  and  temperature  de- 
viations of  a few  degrees  Centigrade,  have  influence  in  changing  long-term 
weather  patterns  in  regions  of  the  globe  removed  from  the  anomaly. 

For  a large  number  of  OTPP's  operating,  one  could  begin  to  examine 
the  effects  of  limited  area  sea-surface  temperature  anomalies  in  various 
locations  downwind  from  the  anomaly.  Since  the  transport  of  atmospheric 
effects  of  the  anomaly  will  depend  ultimately  on  the  large-scale  global 
wind  systems,  we  must  concern  ourselves  with  global  weather  models,  and 
their  ability  to  treat  air-sea  exchange  processes. 

We  will  be  interested  in  four  models  that  have  been  designed  for 
the  prediction  of  global  weather  patterns,  and  have  also  been  used  for 
climate  simulation.  These  GCM's  are  the  models  respectively  of:  Holloway 
and  Manabe  (1971),  developed  at  GFDL  of  NOAA;  Mintz  and  Arakawa  at  UCLA, 
documented  by  Arakawa  et  al  (1969)  and  Gates  et  al  (1971);  Somerville  et 
al. (1974)  at  GISS  (NASA  Goddard  Institute  for  Space  Studies,  New  York); 
and  Kasahara  and  Washington  (1971)  at  NCAR. 
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In  many  ways  these  GCM's  are  quite  arbitrary  in  how  certain  physical 
processes  are  parameterized,  and  the  recipes  vary  with  the  model.  Still, 
these  models  provide  the  best  apparatus  for  performing  controlled  experiments 
on  the  highly  nonlinear  dynamics  of  the  atmosphere. 

For  the  sake  of  brevity,  the  four  GCM's  could  be  partly  character- 
ized with  how  each  treats  items  a)  to  f)  listed  on  page  46.  One  other 
notable  characteristic  in  distinguishing  between  these  models  is  the  choice 
of  Independent  variable  for  the  vertical  coordinate.  The  height  z,  the 
pressure  p,  log  p,  or  the  ratio  a = p/p^  (where  p^  is  the  variable  sur- 
face pressure)  have  all  been  used  by  various  modelers  of  atmosphere  phenomena. 

All  except  the  NCAR  model  use  the  0 coordinate,  with  NCAR  using 
the  height  z.  The  GISS  (9  vertical  layers)  and  the  UCLA  (2  layers)  models 
use  Arakawa's  parameterization  of  cumulus  clouds,  whereas  NCAR  (6  layers) 
and  GFDL  (9  layers)  use  " moist  convective  adjustment".  In  both  the  GISS 
and  UCLA  models  there  is  no  horizontal  diffusion,  whereas  NCAR  and  GFDL 
use  the  nonlinear  viscosity  deduced  by  Smagorinsky  (1963).  The  vertical 
diffusion  is  mixing-length  dependent  in  the  GFDL  model,  while  in  the  GISS 
and  NCAR  models  it  is  given  by  an  expression  (Deardorff,  1967)  containing 
Richardson  number  and  the  potential  temperature  gradient. 

The  Kurlhara  (GFDL)  model  uses  the  box  method  for  generating  the 
global  mesh  (successive  latitude  belts  containing  fewer  boxes)  such  that 
all  boxes  on  the  sphere,  except  at  the  poles,  have  the  same  size.  This 
method  is  a much  more  efficient  way  of  covering  the  globe  than  the  latitudm- 
longltude  intersections,  which  put  too  many  points  at  the  poles.  All  four 


} 


) 


47 


models  use  leap-frog  time  differencing  for  the  advectlon  and  pressure  terms, 
and  a lagged  time  step  for  the  friction  term  (Williams,  1969).  The  time 
step  allowed  by  these  methods  is  very  small,  usually  because  of  the  fast 
gravity  waves  present  in  the  system;  conversion  to  the  semi-implicit  scheme 
of  Robert  (1971)  would  be  highly  desirable. 

Because  of  the  limitations  of  only  two  layers  in  the  vertical, 
the  UCLA  model  does  not  perform  as  well  on  the  baroclinic  aspects  of  the 
problem  as  the  others.  In  overview,  probably  the  GISS  and  GFDL  models 
with  nine  layers  yield  the  best  simulations  of  the  atmospheric  dynamics. 

All  GCM's  are  accompanied  by  a "predictability  limit"  due  to 
imperfect  initial  conditions,  finite  spatial  resolution,  accumulated  time 
truncation  and  perhaps  round-off  errors.  A small  change  in  the  system  of 
equations  or  in  initial  conditions  will  produce  changes  in  the  atmospheric 

j 

variables  which  increase  by  an  order  of  magnitude  every  five  days  or  so;  • i 

thus  one  cannot  readily  distinguish  between  effects  of  small  errors  or  j 

small  changes  in  the  environmental  parameters  (e.g.  solar  radiation,  . 

polar  ice  cover,  etc).  Nevertheless,  since  in  climate  studies  one  is  only 
interested  in  certain  averaged  states  (or  statistical  behavior)  of  the  ^ 

atmosphere,  several  workers  carried  out  GCM  experiments  to  study  the 
seasonal  variations  of  the  atmosphere  or  the  climate  determined  by  dif- 
ferent sets  of  the  environmental  parameters.  The  results  indicate  that 
large-scale  changes  of  the  environmental  conditions  will  indeed  lead  to  a 
computable  climate  change.  A reproduction  of  the  climate  for  one  month  of 
the  year  was  performed  by  Manabe  et  al.(1970),  Washington  & Kasahara  (1970), 
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Gates  (1972),  and  Somerville  et  al.  (1974).  All  chose  the  month  of  January 
and  all  succeeded  In  simulating  some  of  the  essential  features  of  the  zonally 
averaged  mean  flow  field  and  tne  heat  and  water  budgets  (see  figure  13).  In 
all  these  experiments  the  solar  radiation  and  the  sea-surface  temperature 
were  held  constant  at  values  corresponding  to  January. 

The  seasonal  variation  of  the  atmosphere  has  been  simulated  by 
Manabe  et  al  (1974)  who  used  an  eleven-layer  GCM  model  to  study  the  response 
of  the  atmosphere  to  insolation  (incident  solar  radiation)  and  sea-surface 
temperature  over  a three  year  period.  The  model  accurately  simulated  the 
seasonal  variation  of  the  location  of  the  tropical  rain  belt.  Over  the 
oceans  this  belt  had  a tendency  to  form  away  from  the  equator,  which  the 
authors  attributed  to  the  equatorial  belt  of  low  sea-surface  temperatures. 

The  seasonal  variation  of  the  tropical  atmosphere  due  to  seasonal 
changes  in  the  sea-surface  temperatures  has  also  been  investigated  by  Pike 
(1972)  using  a joint  atmosphere-ocean  model.  He  used  a simple  zonally 
averaged  two-dimensional  model  but  Included  two  layers  of  the  upper  ocean 
without  continents.  He  obtained  qualitative  agreement  on  the  migration  of 
a single  mode  of  the  ITCZ  (inter-tropical  convergence  zone)  and  the  lag  of 
the  maximum  subequatorial  sea-surface  temperature  behind  the  overhead  sun. 

A combined  global  ocean-atmosphere  simulation  was  carried  out  by 
Manabe  (1969)  to  Identify  various  influences  of  the  circulation  of  the 
oceans  on  the  general  circulation  and  the  liydrologic  cycle  of  the  atmosphere. 
The  integration  was  performed  for  one  atmospheric  model  year  and  a hundred 
oceanic  years,  the  two  time  scales  overlapping  to  make  the  computational 
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Fig.  17.  Computed  and  observed  January  zonal  mean  fields  of  icni|KTaturc  (‘‘O.  Observed  field  is  l»a»td  on  data 
from  Oon  am)  I\a:jniu.<>on  Negative  regions  are  shaded. 


LATITUDE  (DEG) 


COMPUTED 


LATITUDE  (DEG) 

Fic.  14.  Compuud  ami  «d»sorve<l  January  zonal  mean  fields  of  zonal  wind  (n>  Het  *).  OI»serve<l  field  is  lusid  on  <lata 
from  (Xirt  and  k.i<nius.Min  (ld7l).  \ej»ati\r  (easterly)  regions  arr  »hade«l. 

Fig.  13  - ResulU  of  GISS  general  circulation  model  (GCM)  calculations  of  a January  climate, 
compared  to  observations  [Somerville  et  al.,  1974] 




time  reasonable.  Even  though  the  joint  model  consumed  500  hours  of  CDC 
6600  machine  time,  without  the  deeper  layers  of  the  ocean  coming  into 
equilibrium,  some  of  his  findings  were:  a)  the  reduction  of  the  poleward 
transport  of  latent  and  heat  energy  by  the  atmosphere  due  to  ocean 
currents;  b)  an  increase  in  the  temperature  contrast  between  oceans  and 
continents  in  middle  and  high  latitudes,  enhancing  the  development  of 


cyclones  oft  the  East  Coast;  c)  at  the  equator  the  upwelling  of  cold 
water  tends  to  suppress  tropical  rainfall  over  the  ocean  but  increase  it 
over  the  continents;  d)  the  middle  latitude  rain  belt  shifts  poleward 
due  to  the  ocean  current. 


Next,  we  must  mention  an  interesting  climatic  change  experiment 


by  Washington  (1972),  who  investigated  the  possible  climatic  changes  caused 


by  the  release  of  all  thermal  pollution  associated  with  the  energy  use 


of  all  mankind.  Assuming  a population  of  20  billion  and  a personal  use  of 
15  kw  with  a realistic  population  distribution,  he  found  Lljat  tlie  addition 
or  subtraction  of  3 x 10*“*  kw  to  the  atmosphere  produced  the  same  changes 


in  the  final  state  after  40  days  of  iteration  that  the  addition  of  a small 
initial  error  has  made.  It  appeared  therefore  that  the  effect  of  the 
thermal  pollution  was  smaller  than  the  natural  fluctuations  of  the  model. 
The  result  of  Washington  does  not  preclude  a possible  climate  cliange  by 
thermal  pollution  of  tiie  ocean,  or  at  least  of  its  deep  layers,  because 
the  previous  experiments  of  Pike  (1972),  Manabe  et  al  (1974),  and  Spar 


(1973  a,  b,  c)  have  all  stiown  regional  climate  changes  due  to  local  sea- 
surface  temperature  changes.  Somerville  (1976b)  lias  demonstrated  global 
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changes  in  a study  of  GCM  sensitivity  to  sea-surface  temperature. 

The  GCM's  have  also  been  used  to  study  effects  on  climate 
when  much  larger  changes  were  imposed  through  the  boundary  conditions. 
Williams,  Barry  & Washington  (197A)  have  used  the  NCAR  GCM  to 
simulate  typical  winter  and  summer  weather  conditions  with  a snow 
and  ice  cover  appropriate  to  a glacial  maximum.  Other  studies 
consider  the  effects  of  modifying  the  solar  flux  constant,  as  in 
work  of  Somerville  et  al.  (1976a)  using  the  GISS  GCM  code.  Although 
such  climatic  research  is  not  directly  applicable  to  OTPP  problems, 
it  does  serve  to  illustrate  some  of  the  more  extreme  climate 
modification  studies  based  on  GCM's. 
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6. 


MARINE  ECOSYSTEMS 


The  ecological  balance  processes  in  the  deep  ocean  arc  complicated 
and  have  a large  literature.  Studies  on  the  modeling  of  these  processes 
have  been  given  by  Platt  & Denman  (1975)  and  Raymont  (1966).  Briefly, 
phytoplankton  growth  depends  on  light  intensity  and  the  presence  of  suf- 
ficient nutrient  (nitrate,  ammonia  and  phosphate).  Zooplankton  consume 
phytoplankton,  and  fish  consume  both.  Dead  organisms  tend  to  sink,  and 
then  rot  at  greater  depths,  so  that  nutrient  is  continually  carried  away 
from  the  sunlit  area  by  falling  detritus.  Thus  in  most  of  the  ocean, 
marine  life  is  nutrient-limited.  (Verhoff,  1971;  Thomas  & Dodson,  1968; 
Thomas,  1970;  Macisaac  & Dugdale,  1969;  Corner  & Davies,  1971). 

To  obtain  empirical  equations  describing  these  different  processes, 
together  with  horizontal  and  vertical  advection  and  vertical  turbulent 
transport  of  the  different  organisms,  is  relatively  easy.  To  obtain 
equations  giving  good  agreement  with  observations  over  large  ocean  regions 
is  more  difficult,  and  it  is  fair  to  say  that  only  moderate  progress  has 
been  made  so  far  (Steele,  1959;  Riley,  Storamel  & Bumpus,  1949). 

The  importance  of  upwellings  to  the  marine  ecology  has  long  been 
recognised  (Bogdanov,  Sokolov  & Khroraov,  1968;  Cushing,  1969;  Dugdale, 
1967;  Glooschenco,  Curl  & Small,  1972).  O'Brien  & Wroblewski  (1973), 
Wroblewski,  O'Brien  & Platt  (1974),  and  Wroblewski  & O'Brien  (1975) 
have  modeled  the  dynamics  as  well  as  the  ecology,  in  order  to  explain 
observed  features  of  the  plankton  distribution. 
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In  regard  to  OTPP  operation,  the  limiting  nutrient  is  a passive 
scalar  in  the  near-field,  since  the  time-scale  is  too  short  for  significant 
changes  in  the  marine  life.  But  a detailed  ecosystem  model,  possibly  of 
a horizontal  average  type,  is  required  to  estimate  the  possible  economic 
benefits  produced  by  the  huge  upward  nutrient  flux  associated  with 
OTPP  operation. 
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